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Abstract 

The measurement of Aa/a from the QSO absorption lines of Fe II taken at two different 
\ redshifts (Quast et al, Levshakov et al) is compared with the result based on the 23 data 

' points also obtained from VLT/UVES (Chand et al), to demonstrate that the apparent 

oscillatory time-dependence fitting the latter taken at face value is reinforced by the 



former claimed to be of much better accuracy. 

Results of the two groups for the possible time-variability of the fine-structure constant derived 
from the QSO absorption lines appear to be incompatible with each other. The Keck/HIRES group 
reports a non-zero result, Aa/a = (-0.573 ± 0.113) x 10~ 5 (MFWDPW pQ), whereas one of the 
VLT/UVES measurements produces an apparently null result, (—0.06 ± 0.06) x 10~ 5 for the same 
quantity (case 1 of CSPA 0). Another VLT/UVES group shows (0.04 ± 0.15) x 10" 5 (QRL and 
LCMD 0), thus favoring case 1 of CSPA (to be called CSPA1 in what follows). All of these groups 

O ' 

relied on the many-multiplct method and analyzed the results in terms of the weighted-mean estimate, 
O !~ which is interpreted as a 1-parameter fit in terms of a constant Aa/a. 

We proposed, in contrast, an oscillatory behavior for Aa/a as a function of time based on the 
scalar field supposed to be responsible for the acceleration of the universe 5 ; the required "tracking" 
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behavior can be realized by the scalar field being trapped to a potential minimum, thus causing a 
damped oscillation. The oscillation is supported also by the Oklo constraint indicating a much smaller 
Aa/a at the look-back time s » 0.142 00, somewhat closer to the present time than the range of 
the QSO observation, 0.3 <, s < 0.8. The meteorite constraint at s ~ 0.34, however, is not included 
because of the likely model-dependence 0. From a phenomenological point of view, we tried to fit 
the QSO data in terms of a damped oscillator with 3 parameters. 

We applied the 3-parameter fit in this sense to MFWDPW 0HU] and CSPA [TT] . For the former, 
we found nearly the same degree of fit as in the original weighted-mean analysis, whereas we came 
across a nonzero fit to CSPA1, even with a considerably "improved" Xrd> though its value much smaller 
than 1 might indicate a poor statistical quality of the data used in this particular manner. 

It can be argued 0, on the other hand, that the real accuracy in the analysis in CSPA may not 
have reached the level as presented in 0. It may not be sufficient to detect detailed structure of 
the behavior, if any, of Aa/a. We pointed out JT], on the other hand, that choosing the 1- or 3- 
parameter fit depends on what type of distribution, flat or oscillatory, we expect to be approximated 
from a theoretical point of view. An oscillatory behavior might be easily averaged out by the assumed 
flat 1-parameter fit, for example. In Ref. JT], we in fact demonstrated that the CSPA1 data allows a 
nonzero oscillatory fit not in contradiction with a null result which could be derived by a 1-parameter 
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fit. In this sense, the oscillatory fit may deserve further scrutiny for its own, not being limited strictly 
to the parameters derived currently from CSPAl. 

In this article we apply the same 3-parameter approach to QRL and LCMD, finding that these 
results, claimed to have been analyzed much more carefully than the others, appear to provide a strong 
support for the oscillatory nature of CSPAl. It even seems as if the accuracy of the measurement 
and analysis has been already sufficient to describe the physical content correctly. This might be 
too optimistic a view from current efforts to remove many kinds of systematic uncertainties by using 
larger telescopes, for example, not to mention a very limited number of data points utilized to draw 
our conclusion. We still emphasize that our numerical result obtained from the available data appears 
too impressive to be ignored as a mere coincidence, as long as we take the CSPAl result at face value. 

We start with our 3-parameter fit of a damped-oscillator to y(s) = (Aa/a) x 10 5 as a function of 
the look-back time s = 1 — t/to with to for the present age of the universe [TU] : 

y(s) = a (e 6 ^- 1 ) cos(v - Vl ) - e~ b cos(wi)) , (1) 

where v/s = Vx/s\ — i> kio/sokio = 27rT _1 with the coefficient v± determined by 

vi = tan" 1 ((e-^'o - cos(w oklo )) / sin(w ok i )) , (2) 

due to which y(s) vanishes at s = (today) and s — s kio ~ 0.142 corresponding to the Oklo time 
about 1.95 x 10 9 Gy ago. We modified the old definition a Q id of the amplitude in [jU^HE] by a — a \de b 
for an easier representation of Fig. 2 in the following. 

The upper panel of Fig. 1 for the CSPAl result has been borrowed from Fig. 3 of (TT| with 
a = —0.050 replaced by a x e 31 = —1.11. The redshift z in the observation has been translated to 
s by assuming the cosmological expansion parametrized by to — 1.38 x 10 10 y, h — 0.7, 12a = 0.7. In 
addition to the natural zero at s = 0, we chose the vanishing at s = s kio supposed to be a good 
approximation for much smaller Oklo constraint \y\ < 10~ 2 7 , which likely remains robust even with 
the suggested non-Maxwell-Boltzmann neutron spectrum |12| . We found a considerably "improved" 
result Xrd ~ 0-527 compared with the original 0.95 for the weighted mean in |2J. 

In spite of a critical view against Xrd *C lj we notice that the calculated x?d i s small simply 
because our fitting curve passes nearly through a number of data points particularly in the range 
0.60 < s <, 0.78 corresponding to the range of the redshift 1.2 < z < 2.4, where the data shows a 
pronounced "peak." It is tempting to accept this peak as representing something real, though its 
legitimacy depends crucially on the assumed "sufficient" accuracy reached already, as discussed before. 

The authors of 0] chose Fe II intending to overcome uncertainties coming from the isotopic 
abundances and the radial velocities, among other things. By careful analyses they obtained the 
results for the two redshifts 1.15 and 1.84 with the statistical and systematic errors shown separately. 
They combined, however, the two results to obtain the overall weighted mean y — 0.04 ± 0.15 s t a t, 
which is in favor of CSPAl, but not MFWDPW. From our point of view of the oscillating behavior, 
we attempt to consider their two results separately instead of combining them together. 
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Figure 1: The best 3-parameter fits for CSPAl only (upper panel), CSPA1 + QRL + LCMD (middle 
panel), and CSPAl + LCMD (lower panel), versus the fractional look-back time s as well as the redshift 
z. The curves below s = 0.2 are magnified as indicated in each panel, including the (assumed) zero 
at Sokio = 0.142. The points for QRL, LCMD are marked by crosses with the error bars denoted by 
thicker lines. The horizontal dashed line stands for —0.06 for the weighted mean derived in CSPAl. 

We first consider the 25 data points consisting of CSPAl, QRL( y = —0.04 ± 0.19 s tat at z = 
1.15, s = 0.60) and LCMD(y = 0.24 ± 0.38 sta t at z = 1.84, s = 0.725). Our best 3-parameter fit is 
shown in the middle panel of Fig. 1, together with the corresponding 68%-confidence region, also in 
the middle panel of Fig. 2. The best value of T turns out to be nearly the same as in the upper panel. 
The values of 6, and hence of a are somewhat different. We also find Xrd ~ 0-578 which is larger 
than 0.527 for the pure CSPAl slightly. This suggests that the Fe II data shares some similarity with 
CSPAl. 

We find in fact that the LCMD data point, in particular, lies in the midst of the surrounding 
points that constitute an obvious "peak" of CSPAl. It appears as if the LCMD result obtained with 
much more care than CSPAl reinforced the most important part of the oscillating nature as was 
observed in CSPAl. In order to support this view further, we applied the same 3-parameter fit to 
the 24 data points including the CSPAl and LCMD data, as shown in the lower panels of Figs. 1 
and 2. The similarity with the pure CSPAl is even more obvious, resulting in the still smaller value 
Xr d = 0.522. This implies, on the other hand, that the QRL result contributes a little to make the 
agreement poorer. We still point out that the fitting curve in the middle panel of Fig. 1 passes barely 
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Figure 2: The 68% confidential volumes in 3-dimcnsional a,6, T space are shown as 2-dimensional 
cross sections for each value of b, as indicated alongside the corresponding contours. The contours 
drawn by solid curves also with the crosses inside correspond to the best values of b, whereas those 
by dotted curves with filled circles are for representatives of different values of b. The three panels 
correspond to those in Fig. 1, respectively. In the middle panel, the horizontal scale has been stretched 
by the ratio 0.3 and shifted by 0.2 to the right. 

off the 1 standard deviation at z = 1.15 (QRL), implying the "damage" to be by no means serious. 

We also add that the result of CSPA1 together with QRL alone can be fitted by a = —2.18, b = 
5.9, T = 0.136, Xid = 0.595. As anticipated, the reduced chi-squared is a little larger than that 
corresponding to the middle panel of Fig. 1, but not significantly again. 

To summarize, we emphasize that the Fe II samples due to QRL and LCMD, taken at face value at 
the separate z values, reinforce the oscillatory behavior as illustrated in CSPA1. We admit, however, 
that the strongest support comes from a single data point at z — 1.84 (LCMD). We still notice that 
this data point happens to occur at z around which we find the peak behavior crucially important to 
reveal the oscillatory structure in CSPA1. This latter analysis might be subject, however, to future 
improvement with respect to the accuracy, as was argued in the most recent version (v3) of LCMD. 
We must be open to a more conventional wisdom that the required accuracy has not been fully realized 
yet to establish the presence of a peak unambiguously. In order to determine if the oscillatory time- 
depcndcncc of a is in fact justified, we urge that an emphasis should be placed in the future analysis 
to make careful measurement at different z values even for relatively sparse data points. 

We close the paper by adding a theoretical remark. Once we accept an oscillating behavior, it is 
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irrelevant to ask if a was larger or smaller in the past. As a consequence, rather than dividing the 
past flat value of y obtained from the weighted-mean by an averaged time span, we estimate today's 
rate of change of a in terms of the tangential slope of y(s) at s = 0; (a/a)o = —t^ 1 y'(0) x 10 -5 
The upper and the lower panels in Fig. 1 yield a particularly "large" value, —0.96 x 10 5 y _1 , which 
is close to the upper bound ~ 2 x 10~ 15 y _1 reached recently by the laboratory experiment |13|. 

I would like to thank Paolo Molaro for his valuable remarks on the nature of the data of QRL and 
LCMD. I also thank Aldo Fiorenzano and Naoto Kobayashi for their illuminating discussions on the 
subject. 

References 

[1] M.T. Murphy, V.V. Flambaum, J.K. Webb, V.V. Dzuba, J.X. Prochaska, and A.M. Wolfe, Proc. 
Astrophysics, Clocks and Fundamental Constants, 16-18 June 2003, Bad Honnef, Lect. Notes in 
Phys. 648, 131 (2004), astro-ph/03i03T8| 

[2] H. Chand, R. Srianand, P. Petitjean, and B. Aracil, Astron. Astrophys. 417, 853 (2004): R. 
Srianand, H. Chand, P. Petitjean, and B. Aracil, Phys. Rev. Lett. 92, 121302 (2004). 

[3] R. Quast, D. Reimers and S.A. Levshakov, A & A, 415, L7 (2004). 

[4] S.A. Levshakov, M. Centurion, P. Molaro and S. D'Odorico, astro-ph/0408188. 

[5] Y. Fujii and K. Maeda, The scalar-tensor theory of gravitation, Cambridge University Press, 2003. 



[6] A.I. Shlyakhter, Nature 264, 340 (1976); unpublished (1983), physics/0307023, T. Damour and 
F.J. Dyson, Nucl. Phys. B480, 37 (1996). 

[7] Y. Fujii, A. Iwamoto, T. Fukahori, T. Ohnuki, M. Nakagawa, H. Hidaka, Y. Oura and P. 
Moller, Nucl. Phys. B573, 377 (2000); Proc. ND2001, 7-12 October 2001, Tsukuba, Japan, 



hep-ph/0205206 



[8] Y. Fujii and A. Iwamoto, Phys. Rev. Lett. 91, 261101 (2003). 
[9] Y. Fujii, Phys. Lett, B573, 39 (2003). 

[10] Y. Fujii, Proc, Astrophysics, Clocks and Fundamental Constants, 16-18 June 2003, Bad Honnef, 
Lecture Notes in Physics, Springer; hep-oh/0311026 



[11] Y. Fujii and S. Mizuno, astro-ph/0404222 



[12] S.K. Lamoreaux and J.R. Togerson, Phys. Rev. D69, 121701 (2004). 

[13] E. Peik, B. Lipphardt, H. Schnatz, T. Schneider, Chr. Tamm and S.G. Karshenboim, Phys. Rev. 
Lett. 93, 170801 (2004). 



5 



